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Abstract 

This paper is devoted to the transport process of momentum, the 
viscosity. We outline the historical development of hydrodynamics, 
the branch of physics in which the shear and bulk viscosities are ex- 
pressed. We show the transition between the mechanical philosophy 
claiming that the conservation of energy is coupled to cyclic, reversible 
dynamics, and the birth of a new science, thermodynamics, in which 
the conservation of energy (the first law) is associated with the irre- 
versibility (the second law). Through the 19th century, the kinetic 
theory of gases was developed and accumulated successes establish- 
ing relationships between the microscopic dynamics of the atoms and 
molecules, and macroscopic properties, among which the viscosity. In 
this context appeared the paradox of irreversibility to which Boltzmann 
was the first to provide an explanation based on a statistical interpre- 
tation of the entropy. We finally outline the history of the theory 
of chaos and show how the chaotic character of the microscopic dy- 
namics provides a dynamical explanation of the irreversible processes 
such as transport, and in particular, the viscosity. In this context, we 
present three successfull approaches developed the last two decades: 
the thermo stated- system method, the escape-rate formalism and the 
hydrodynamic-mode method. 

*This paper consists in a revision of the introduction of the PhD thesis presented 
on the 21th of September, 2005, entitled : Viscosity and Microscopic Chaos : The 
Helfand-moment Approach (Viscar dy, 2005| > . 
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1 Introduction 



The nineteenth century was rich in scientific discoveries and in new con- 
cepts. In particular, as regards the subject discussed in this paper, the nine- 
teenth century appears as a turning point. Indeed, the industrial revolution 
gave rise to the development of a new science, thermodynamics, which em- 
phasized the irreversible character of macroscopic phenomena (such as heat 
dissipation). On the other hand, the development of chemistry and, later, of 
the kinetic theory of gases, transformed the concept of matter from a contin- 
uous to a discrete description. In other words, matter should be considered 
as made of very small particles: the atoms and the molecules. However, the 
motion of these entities is governed by the reversible equations of Newton. 
As a consequence appeared the famous paradox opposing macroscopic and 
microscopic descriptions, which has given rise after the seminal works of 
Maxwell and Boltzmann to a large number of contributions addressing this 
paradox. 

Maxwell was the first to introduce viscosity in the framework of kinetic 
theory and nonequilibrium statistical physics. Since then, it has remained 
a driving concept in these fields. We therefore devote the paper to the case 
of viscosity as a transport property. 

In the sixties a new revolutionary field of the science of the twentieth 
century was born : Chaos. Sometimes compared to the scientific revolu- 
tions implied by quantum mechanics and the relativity, the phenomenon of 
chaos has been discovered in most of the natural sciences such as physics, 
chemistry, meteorology and geophysics. Our intent is to show the role of 
microscopic chaos in statistical mechanics of irreversible processes. 

The paper is organized as follows: we will start in section 2 with the 
presentation of the discovery of the concept of irreversibility as a fundamen- 
tal property in nature. The development of hydrodynamics as the domain 
of physics being particularly concerned by the process of viscosity will be 
considered in section 3. Further, in section 4, we focus on the kinetic theory 
that has a long history and has played an important role in the establish- 
ment of the statistical mechanics, that is to say, in the understanding of 
the relation between the microscopic and macroscopic behaviors. Finally, 
in section 5, we cover the brief history of the development of the theory of 
chaos and of dynamical systems which provides a better understanding of 
the origin of the irreversible processes. 
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2 Irreversibility 



The history of modern science holds its richness in the variety of concep- 
tions which were developed successively throughout the years 1 . In the 17th 
century, it was the doctrine of the clockmaker God or mechanical philosophy 
which prevailed above the others. 

The development of these ideas began with the work of Boyle (1627- 
1691) and others. According to his doctrine, Boyle stated that the world 
works like a clock: once created and energized, it can run forever in a de- 
terministic way and without any need for a divine intervention. To ensure 
that this "clockwork universe" never runs down, Descartes (1596-1650) more 
or less intuitively introduced the statement that the total amount of mo- 
tion in the world must remain constant. He defined this quantity as the 
scalar momentum mv (quantity of matter multiplied by its scalar speed 2 ) 
(Blackwell, 1966). However, as Descartes himself later observed, experi- 
ments did not confirm his enunciation of the conservation law of motion. 
Huygens (1629-1695) corrected it in 1668 by modifying it into the vector 
form. Huygens worked also on the problem of collisions. This revised law 
allowed him to claim that the vector sum of the product of the mass by 
the vector velocity mv remains unchanged after a collision even if it was 
inelastic and with dissipation of energy after the collision itself. Introducing 
this modification, it then appears that the world could stop after a certain 
amount of time. This was totally contradictory with the clockworld concept. 

The only way to avoid this possibility was to postulate that the matter 
is composed of elastic particles. If the macroscopic objects lose some motion 
after a collision, it is only in appearance because the motion is transferred 
to the invisible particles of the objects. Therefore, the idea emerged that 
heat is related to a rapid motion of the invisible parts inside the macroscopic 
bodies. 

In the 18th century, the study of the motion of bodies and its change 
after collisions was of great importance. Moreover, in addition to this study, 
another quantity so-called vis viva, mv 2 , was introduced for the first time. 
Besides, its conservation law in elastic collisions was remarkably proved by 
the same man who helped formulate the conservation of momentum: Huy- 
gens. This major contribution induced the first step to the central quantity 

1 For general informations on this section, see for example Holton and Brush (1985) 
and Brush (1983) . 

2 More precisely, Descartes considered the volume of the bodies instead of their mass, 
the volume being for him the real measure of quantity of matter to be considered, as 
pointed out by Blackwell (1966) . 
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acknowledged in mechanics: the energy. The concept of vis viva was used 
up to the 19th century, when the factor ^ was added and it then became 
"kinetic energy". 

The world-machine concept is often attributed to Newton (1642-1727) 
because of the importance played by his brilliant Principia published in 1687 



(Newton, 1999 ). In his masterpiece, the Newtonian dynamics is depicted for 



the first time although it clearly appeared that the author was really opposed 



to it as he pointed out in his Opticks (Newton, 1952 1 published in 1704 



By reason of the [...] weakness of elasticity in solids, motion is much more 
apt to be lost than got, and is always upon the decay. 

Therefore, in a certain sense, he suggested the existence of the dissipation of 
motion in the world. However, it should be specified that this viewpoint has 
to be considered in the light of the polemic between Newton 3 and Leibniz 
concerning the role of God (Koyre, 1957). 



During the 18th century the concept of the "Newtonian world-machine" 
continued to dominate the thought of scientists. This century is also the 
one of the birth of Geology. It began with the question of the temperature 



of Earth's interior ( Brush, 1994 1. Yet in 1693 Leibniz thought that Earth 
was originally hotter and cooled down, at least on the outside. Later Buffon 
(1707-1788) studied this question by leading some experiments on the cool- 
ing of heated spheres of iron. By considering a molten Earth at the origin 
he hence found that our planet was about 75,000 years f de Buffon, 1774 ). 



At the opposite, the founder of Geology James Hutton (1726-1797) dis- 
agreed with this theory of cooling. Defending his "Uniformitarism" (princi- 
ple telling that the geological processes in the past have to be explained by 
using only the laws and physical processes that can now be observed) and 
accepting the hypothesis that the interior of the Earth is much hotter than 
its surface, he thought this situation had been like that forever. For him the 
geological processes are cyclic: alternance of periods of erosion and denuda- 
tion implying the destruction of the mountains, and periods of uplifts of new 
continents (thanks to the subterranean fires). Actually one of his disciples 
Playfair promoted the uniformitarist viewpoint by citing the mathematical 
works by Lagrange and Laplace showing the cyclic movement of the planets 
around the sun. This position was perpetuated by Lyell (1797-1875) during 
the 19th century. 

In the early 19th century, the first mathematical theory describing the 
propagation of heat was developed. Indeed, the French mathematician 

3 More precisely, his student and faithful friend Clarke (1675-1729). 
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Fourier (1768-1830) was interested by the problem of the cooling of the 
Earth. This was the main motivation that lead him to study the heat con- 
duction in solids. In 1819, in his Memoire sur le refroidissement seculaire du 
globe terrestre (Fourier, 1819), he came up with an equation which presented 
a significant feature: unlike Newton's laws of motion, Fourier's equation is 
irreversible. Any system in which a temperature difference exists presents 
a heat flow from the high to low temperature. This was the discovery of an 
equation describing processes in which the past does not play the same role 
as the future. His theory therefore presents an important turning point in 
the history of physics due to the powerful analysis he developed (what we 
now call Fourier's analysis), and because it is explicitly based on a postulate 
of irreversibility. In the first half of the 19th century the eventual contradic- 
tion between the reversible Newtonian dynamics and the irreversible heat 
conduction did not appear because the most studied problem at this time 
was at the phenomenological level and the Newtonian mechanics had already 
been successfully applied to problems with dissipative forces (friction, etc.). 
The contradiction appeared during the first attempt to explain the macro- 
scopic processes in terms of the assumed reversible Newtonian dynamics of 
the particles composing the system. 

A few years after Fourier's work on the propagation of heat, Carnot 
(1796-1832) published in 1824 his essay Reflexions sur la force motrice du feu 
dCarnot, 1824| ). Driven by engineering ideas, he focused his research on the 
question of the limited efficiency of the steam engine. By using the caloric 
theory^ appearing during the 18th century and developed by Lavoisier (1743- 
1794), Carnot concluded that a difference of temperature existing between 
two bodies gives the possibility of doing work by allowing heat to expand 
a gas as the heat flows from the hot body to the cold one 5 . However, 
an engine wrongly designed results in loosing the chance of doing work that 
might have been done. Therefore, proving that heat always flows from hot to 
cold bodies in the engine as it does anywhere in nature, and discovering that 



4 The caloric theory supposed that heat is a fluid composed of particles independent 
of the rest of matter. These particles repel each other but are attracted to the particles 
of ordinary matter. According to this theory it has to be an eventual material substance 
and thus should be limited. Hence the increase of the temperature of water by rotating a 
drill should be due to the transfer of the so-called caloric from the drill to the water. But, 
in 1798, Thompson (1753-1814) realized that one can produce heat without any limit. He 
then concluded that heat is not a chemical substance or a material substance but is the 
expression of a movement. However, since he did not propose any alternative theory of 
heat, the caloric theory had an important influence until the 1830's and the works done 
in particular by Joule. 

Actually later, Carnot abandoned the caloric theory. 
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this equalization represents a loss of the opportunity to produce mechanical 
work, Carnot pointed out the existence of dissipation. More important, the 
discovery made first by Carnot in 1831, then Mayer (1814-1878) in 1842 
( Mayer, 18421 and later by Joule (1818-1889) ( [Joule, 18471 ), brought up a 
fundamental statement: work and heat are actually two different expressions 
of the same quantity: the energy. Indeed, although heat and mechanical 
work seem to be two independent quantities, we now understand that heat 
can be transformed into work and vice versa. In other terms, the energy of 
an isolated system remains constant. A science of transformations was born: 
thermodynamics. This equivalence between heat Q and work W consists in 
the so-called first law of thermodynamics: the conservation of the internal 
energy E of a closed system: dE = dQ + dW . 

As we have seen, the original problem giving birth to thermodynamics 
is the separation between the concepts of conservation and reversibility. In 
mechanical transformations, the conservation of energy (at the beginning 
the vis viva) coincides with the idea of reversibility. On the other hand, the 
physico-chemical transformations can preserve the energy while they cannot 
be reversed. Hence it became necessary to define a quantity characterizing 
this irreversibility. In 1854 Clausius (1822-1888) introduced the quantity 
dQ/T as a measure of the quantity of work lost during the transfer of heat 
from a hot to a cold body ( Clausius, 1854} |Clausius, 1856 ). Later, in 1865, 
he gave it the famous name entropy 6 ( |Clausiu~ T865J. Thanks to this new 
quantity, the second law of thermodynamics can now be stated very sim- 
ply: the entropy of an isolated system always tends to increase. We thus 
have a well-established principle regarding the irreversible processes. We 
may therefore state that the two laws of thermodynamics combine the two 
concepts of conservation of energy and of nonreversibility of macroscopic 
phenomena. 

As we will later see in section @J the atomic theory increased its influence 
during the 19th century. The development of the kinetic theory of gases 
played an important role to prove the discrete character of matter. One 
of the most important leaders in this area was Boltzmann (1844-1906). He 
particularly proposed a more general definition of entropy, in terms of the 
probabilities of molecular arrangements, that can change even when there 
is no heat flow (e.g. the system become more disordered). His theory based 
on a statistical description was a first understanding of the link between the 

6 The ethymology of this word expressed in the most clearly way the essence of ther- 
modynamics: indeed it comes from the Greek words evep^ua (energy) and Tpoirq (trans- 
formation) . 
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atomic level and the macroscopic phenomena, that is between the reversible 
microscopic dynamics and the macroscopic irreversible processes. 

In the nineteenth-century society, in which social and economical activi- 
ties, science and technology induced a progress never observed before in the 
history of Humanity, arose the idea of the evolution of species with Lamark 
(1744-1829) and later with Darwin (1809-1882). According to Darwin the 
species that we know today are the result of a long process of natural selec- 
tion. This mechanism has induced the production of increasingly complex 
and organized beings, from the unicellular bacteries to the mammals such as 
the human beings. Whereas thermodynamics introduced a quantity measur- 
ing the continuous growth of disorder, biology put in evidence the continuous 
growth of order, of organization in the biological world. As a result, these 
two points of view influenced some philosophers: on the one hand, Adams 
(1838-1918), influenced by the discovery of constant dissipation of energy in 
the universe, expressed a pessimistic vision concerning the future of the soci- 
ety, as argued in The Degradation of the Democratic Dogma ( Adams, "2 004), 
and made the second law of thermodynamics an explicit basis for the ten- 
dency of history (Burich, 1987); on the other hand, the Darwinism strongly 
influenced philosophers like Spencer (1819-1903) who thought that all in uni- 
verse goes gradually from a state of confused homogeneity to a differentiated 
heterogeneous state ( Freeman, 1974} Rousseau, 1945 ). These two points 
seemed to be contradictory. One had yet to wait for the development of the 
nonequilibrium thermodynamics, especially attributed to Prigogine (1917- 
2003) and his coworkers, to resolve this apparent contradiction. Indeed 
they showed that, out of thermodynamic equilibrium, i.e. in open systems, 
matter is able to exhibit self-organization, although the internal production 
of entropy is absorbed by the environment (Glansdorff and Prigogine, T97T1 
Nicolis and Prigogine, 1977). These recent works have pointed out the im- 



portance of nonequilibrium phenomena and have led to important advances 
on the properties of nonequilibrium systems. It is in this general context 
that this thesis is written. 



3 Hydrodynamics and viscosity 

The vital importance of water throughout all civilizations induced a 
special interest for the study of the properties and behaviors of fluids (in 
particular of water). The earliest quantitative application of "real fluid" or 
viscosity effects was by the ancient Egyptian Amenemhet (~1600 BC). He 
made a 7 degree correction to the drainage angle of a water clock to account 
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for the viscosity change of water with temperature (which can be significant 
between day and night in this region). And in general the early centres 
of civilization in Egypt, Mesopotamia, India and China systematically used 
various machines for irrigation and water supply (Blair, 1944 ). In particular, 



Archimedes (287-212 BC), is not only considered as the father of hydrostatics 
for his famous law, but has also developed the so-called Archimedes screw, a 
water elevating machine, which has been used for different purposes (e.g. to 
extract water from rivers). Notice also that Heron of Alexandria (~ 10-75) 
treated in Mechanica the problem of friction between a body and a plane 
and, considering a horizontal, frictionless plane, he said 

We demonstrate that a weight in this situation can be moved by a force less 
than any given force 7 . 

The earliest extensive treatise devoted to practical hydraulics is a man- 
ual on urban water supply written by Roman soldier and engineer Sextus 
Julius Frontius (first century BC), inspector of the aqueducts and the public 
fontains in Rome. He noted that the amount of water discharged through 
an orifice in a given interval of time depends, not only on the size of the 
orifice, but also on its depth h below the surface of the water in the reservoir 
from which the water is supplied. Starting upon the results of Galileo's ex- 
periments with falling bodies, Torricelli (1608-1647) came to the conclusion 
that the velocity v of the water exiting out of the reservoir is proportional to 
the square root of the h. This is Torricelli 's theorem which has been math- 
ematically expressed as v = \J2gh later by Daniel Bernoulli (1700-1782). 
On the other hand, the centuries since the Renaissance have been charac- 
terized by the major influence of the engineerings on the development of 
physics. More specifically, the construction of bridges and canals induced 
a lot of theoretical studies on the flow of fluids. By observing the behav- 
ior of the flow of water in rivers, da Vinci (1425-1519) already came to the 
conclusion that, when the river becomes shallower or narrower, the water 
flows faster. Later Castelli (1577-1644) confirmed this result and gave the 
so-called velocity- area law: vA = const, where v is the velocity of water 
and A the cross-sectional area of the flow. It consists of the first idea of 
continuity of flow (for an incompressible fluid), which was later developed 
by Euler. 

The term hydrodynamics was first used by Daniel Bernoulli as the title 
of his book hydrodynamica published in 1738 (Bernoulli, 17381. His the- 



ory was original because he was the first scientist to combine the sciences 
7 Cited by Russo (2004). 
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of hydrostatics (considering the pressure) and hydraulics (considering the 
motion of fluids). From the conservation of vis viva he reached the famous 
so-called Bernoulli's principle relating the velocity of flow at a point in a 

2 

pipe to the pressure there: ^ — \~ P = const. But in the modern sense, hy- 
drodynamics began with the work of d'Alembert (1717-1783) and especially 
Euler (1707-1783). 

During his Berlin period as an engineer, Euler was in charge of the 
construction of canals, the supply of water for King's Sans-Soucis Palace, 
and the improvement the water turbine. However, his strong interest for 
mathematics, mechanics and physics was much greater. Hydrodynamics was 
among his numerous topics of interest. In 1755, he derived the fundamental 
equations of hydrodynamics (Euler, 1755) by introducing the new and main 
concept of fluid particle 8 . First he proposed the modern and general form for 
the equation of continuity (^Q) , which expresses the conservation of matter 9 

dp dvj , . 

where p is the mass density and Vj the i-component of the velocity of the 
considered fluid particle. Furthermore, considering a small parallelepiped 
and pressure P acting on its different faces, by starting from the Newton's 
second law, he derived the general equations for the motion of an ideal fluid, 
the so-called Euler 's equations of motion 

dvi f dvi dvA dP 

( "W = p ['dt +Vi ar J ) = -aF i +fit (2) 

where fa is the i-component of the external force (gravity,...). 

But these equations do not present the property of fluids which was 
more and more of interest at the beginning of the nineteenth century, that 
is the internal friction, or the viscosity 10 in modern terminology. Indeed, the 
dissipative feature of the fluid dynamics does not appear in Euler's equations 
and it was not until then considered as a fundamental property in nature. 



A fluid particle is imagined as an infinitesimal body, small enough to be treated 
mathematically as a point, but large enough to possess such physical properties as volume, 
mass, density, and so on. 

9 We use the Einstein's convention of summation over repeated indices: «j§^ has to 

be understood as V* . v,- 

10 The word viscosity derives from the latin word viscum for mistletoe. The mistletoe 
berries produce a viscous glue which was used to lime birds. Viscosity is expressed in Pa.s 
or in Poise (10 Pa.s) in the honor of Poiseuille (1797-1869). 
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The first studies on this property go back already to Newton's Book II of 
his Principia ( Newton, 19990 . His theory stated that, if a portion of fluid 
is kept moving, the motion gradually communicates itself to the rest of the 
fluid. He ascribed to this effect the term defectus lubricitatis, that is, a lack 
of slipperiness. It is nothing but viscosity. Considering two particles of a 
viscous flow sliding one over another, then there exists friction, or viscous 
resistance along the surface of contact. The force of this resistance per unit 
area is nowadays known as the shear stress r. The hypothesis proposed by 
Newton 11 is that the shear stress depends on the speed with which the two 
layers slide one over another. By extension, we obtain that the shear stress 
is proportional to the velocity gradient in the perpendicular direction of the 
movement of the fluid layers: r ~ 

At the end of the 18th century and the beginning of the 19th century, 
numerous and extensive investigations on the flow of pipes and open chan- 
nels started to be carried out by scientists and engineers such as Du Buat, 
Girard, Prony, Coulomb, and especially Poiseuille (1797-1869). Being in- 
terested in the circulation of blood in capillary vessels, he then used glass 
capillaries of very much smaller bore that any of his predecessors. Long 
series of measurements of the quantity Q of liquid discharged in unit time 
were carried out in function of the different factors. He then obtained in 
the forties the relation Q = where P is the pressure, D and L re- 

spectively the diameter and the length of the capillary ( Poiseuille, 1 840a; 
Poiseuille, 1840T1 |Poiseuille, 18461 1 . And in 1856, Wiedemann (1826-1899) 



as well as Hagenbach in 1860, deduced mathematically, by making use of 
Newton's hypothesis, that the constant K is such that dWiedemann, 18561 



Hagenbach, 18 60) 

where R is the radius of the capillary. A new factor rj appears which Wiede- 
mann proposed to call Zdhigkeitconstante der Flilssigkeiten, that is, the vis- 
cosity coefficient of the liquids (Wiedemann, 1856). So for the first time the 
term viscosity was used in a technical sense. 

In parallel, the first scientist having theoretically used the property of 
viscosity in the fundamental equations of hydrodynamics was the French 

11 In his Principia, Newton made precisely the following hypothesis: 

The lack of slipperiness of the parts of a fluid is, other things being equal, 
proportional to the velocity with which the parts of the fluid are separated 
from one another (Book II, Section 9 in Newton (1999) ). 
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civil engineer Navier (1785-1836) in a memoire read in 1822 to the " Academie 
des Sciences" ( JNavier, 1823| ). His idea was that any pressure tends to reduce 



the distances between the molecules of the fluid. Taking into account the 
intermolecular repulsive forces produced by this action, which gave him the 
explanation of the Newton's lack of slipperiness, he added to the pressure 
acting on a small volume the contribution of this effect (Pdij — ► P5ij — n^), 
that is a dissipative term implying the presence of the viscosity coefficient n. 
Being rediscovered by Stokes (1819-1903) in 1845 without introducing any 
molecular behavior (Stokes, 1845), the history has given the name of these 



two persons to the famous Navier-Stokes equations 

dvj _ dP | / d 2 Vj \ | 
dt dri \ drj drj J 

On the other hand, in 1829, Poisson (1781-1836) derived an equation in a 
memoire devoted to the motion of solid bodies and fluids in which appeared 



a second constant (Poisson, 1831 ), in addition to the viscosity introduced by 
Navier. These two constants measure the response of the substance to two 
different types of forces: the first tending to shear it, and the second tending 
to expand or contract it. We saw above that the first type is characterized 
by the viscosity coefficient r] which is more precisely called shear viscosity. 
But until sixty years ago, most of the standard works did not consider the 
other constant rj' called the second viscosity coefficient by making it simply 
proportional to the shear viscosity. Since Stokes (1845) , it was usually 
admitted that 

n' + -n = o, (5) 

which is generally called Stokes' relation. Stokes used this relation by argu- 
ing that in most cases of practical interest the rate of expansion would be 
so small that the value assigned to this sum would be unimportant. And 
before sixty years ago no direct determination of rj has yet been made. 
Indeed it is only in 1942 that Tisza made this determination by evalu- 
ating the ratio when ultrasonic vibrations are passed through a liquid 



(Tisza, 1942). He then showed that this vanishing relation is not justified 
in general and that the re-introducion of the second viscosity coefficient is 
necessary in order to get an agreement between the theory and experiments 



(Markha m et al., 1951} Karim and Rosenhead, 1952} Brush, 1962b ). Con 



sequently, a new coefficient ( can be defined as 



C = V' + 3»7 (6) 
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and is the so-called bulk viscosity coefficient. 

The bulk viscosity is more difficult to measure experimentally than the 
shear viscosity. Moreover, its measure is generally less known. The different 
ways to measure it imply the phenomenon of absorption of sound in liquids 
and gases. Before Tisza (1942) , it was well-known that the absorption of 
sound in polyatomic gases and in liquids obtained experimentally exhibited 
quite great disagreements with the absorption predicted by the classical 
theory based on viscosity. Previously, it was usually admitted that the 
absorption is given as 



2,, 2 
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2ir*v 



1 

Cy 



n 4 ■ 

c P J 3 



(7) 



where v is the acoustic frequency, p the mass density, v s the acoustic velocity, 
cy and cp the specific heats respectively at constant volume and constant 
pressure, and k the thermal conductivity. Hence, only the shear viscosity 
was considered. In his paper flTisza, 1942D , he pointed out that Stokes' as- 
sumption (JSJ) is not justified except in dilute monoatomic gases. However, 
for polyatomic gases and liquids, the bulk viscosity should be quantitatively 
quite important, sometimes more than the shear viscosity (see Table 
Consequently, Tisza obtained a modified expression for the absorption coef- 
ficient in which the bulk viscosity is added in the expression Q as follows 
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(8) 



or by considering the ratio of the viscosity coefficients ^ 



C _ 4 q obs - q el 
V 3 Q! c i 



4 \cy cp J rj 



0) 



On the other hand, the tangential force acting for example on the plane 
xy should be related to the sliding in the two perpendicular directions com- 
posing this plane, that is + in the direction of x, and -qj + ^f in the 
direction of y. This is why Barre de Saint- Venant (1797-1886) combined both 
of them in 1843 ( de Saint- Venant, 1843 ) and obtained more precise expres- 
sions for the pressure tensor, e.g. P xx = P+2r]^- and P xy = r\ {^f- + ^ J ■ 
It allows one to write the Navier-Stokes equations in a more complete form 
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11 x 1(1^ Pa <s 
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1 Q& 


n 
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Ar 


2.27 





H 2 


0.887 


32 


N 2 


1.66 


0.8 


o 2 


2.07 


0.4 


co 2 


1.50 


1.000 


Air 


1.85 


0.6 


H 2 (liquid) 


85.7 


3.1 


Ethyl alcohol 


110 


4.5 


Glycerine 


134.000 


0.4 



Table 1: Shear and bulk viscosities obtained experimentally at 300 K and 1 
Atm in monoatomic, diatomic and polyatomic gases, and in liquids; values 
given by Thompson (1972). 

as 12 



dvi d 



5-._ ( dvi | dv i 2 § dvi \ c 5 9vi 



+ /,. (10) 



When the fluid is considered incompressible = 0), one recovers the 
original Navier-Stokes equation (|4"|). Let us also mention that these equations 
are valid in so far as Newton's assumption is reasonable. A wide range 
of fluids (water, air, etc.) follows the linear relation between the shear 
stress r and the velocity gradient. But liquids composed by lengthened 
molecules (polymers, etc.) obey to different viscosity laws, which are called 
non-Newtonian. 

As we will see in the next section, the 19th century and the beginning 
of the 20th century before the first world war has seen the decisive develop- 
ment of the scientific atomistic conception of matter and the will to derive 
macroscopic phenomena like viscosity in terms of the Newtonian microscopic 
dynamics. These works gave rise to the powerful kinetic theory of gases. 



12 For further information on the history of hydrodynamics, see Rouse and Ince (1957), 
Tokaty (1994), Mikhailov (1994) and Stokes (1846). A more rigorous derivations of Navier- 
Stokes equations can be found for example in ( |Landau and Lifshitz, 1959| . 
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4 Kinetic theory of gases 



The idea of atomism goes back to the fifth century BC with the ancient 
Greek philosophers Leucippus and Democritus. However, the history of the 
kinetic theory of gases does not really begin until the seventeenth century 
when Torricelli, Pascal (1623-1662) and Boyle first established the physi- 
cal nature of the "sea" of air that exerts mechanical pressure on surfaces 
in contact with it. This development of the concept of air pressure can be 
considered as part of the change in scientific attitudes which led to the expla- 
nations based simply on matter and motion instead of postulating "occult 
forces" or teleological principles. Consequently, Boyle (1662) discovered the 
law PV = const now called Boyle law. 

Usually we consider Daniel Bernoulli as the first scientist having pro- 
posed a kinetic theory of gases. Actually, in his famous treatise Hydrody- 
namica flBernoulli, 17381 ), he gave a derivation of the gas law of Boyle and 
Mariotte from "billiard ball" model, assuming that the gas consists of a very 
large number of small particles in rapid motion. Moreover, using the prin- 
ciple of conservation of vis viva he concluded that pressure is proportional 
to the square of the velocities of gas particles, and thus proportional to the 
temperature. In other terms heat is nothing but atomic motion. However, 
in spite of these results, Bernoulli's model was forgotten. He was about a 
century ahead of his time with his theory. The latter success of Lavoisier's 
caloric theory buried it by proposing another conception of the matter based 
on the equilibrium between the caloric repulsion (due to the atmosphere of 
caloric whose density increases with temperature) and the gravitational at- 
traction between matter particles. 

In chemistry, the atomic theory emerged at the beginning of the 19th 
century with Dalton (1766-1844). Thanks to the laws of chemical discon- 
tinuities, that is the law of constant proportions 13 discovered in 1806 by 
Proust (1754-1822) and the law of multiple proportions 14 in 1802-1804 by 
Dalton himself, he proposed in 1808 a discrete view of matter, which should 
be composed by indivisible entities: the atoms which are characterized by 
different atomic weights. But these laws and experiences at that time did 
not reject a concept of equivalents first proposed in 1792 by Richter (1762- 
1807). Because of different ambiguities amongst the partisans of the atomic 
theory, the equivalentism had a great success during the 19th century. How- 
ever, many new phenomena and laws were discovered that only the atomic 



13 The proportion between two combining elements cannot vary continuously. 

14 When two elements can combine according to different ratios, the ratios of their weight 
in the different cases are simple, that is integer numbers. 
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hypothesis could explain. Among these, the law of Gay-Lussac (1778-1850) 
discovered in 1809 15 , which allowed to Avogadro (1776-1850) to emit his 
famous hypothesis (1811) saying that all gases considered in the same con- 
ditions of temperature and pressure contain, for equal volumes, the same 
number of molecules 16 . We can also cite the law of Dulong (1785-1838) and 
Petit (1791-1820) 17 (1819); the emergence of structural chemistry, in partic- 
ular the development of the theory of valence by Kekule (1829-1896), organic 
chemistry and the concept of isomerism, and stereochemistry. But one of the 
most important works having played a fundamental role in the development 
of the atomic theory and its acceptance is the establishment of the periodic 
table of the elements in 1869 by Mendeleev (1834-1907) flMendeleev, 1869a| 
Mendeleev, 1869b). Indeed, by assuming his periodic law in the properties 



such as the valence of elements when these are listed in order of atomic 
weight, he predicted the existence of unknown elements corresponding to 
the gaps in his famous table and to which a given atomic weight is associ- 
ated (Brush, 1996). The discovery of some of them in the following years 
confirmed the periodic property of elements. In the polemic between atom- 
ists and equivalentists, Mendeleev's brilliant idea heavily contributed to the 
acceptance of the atomic conception of the matter amongst the chemists 18 . 

While the caloric theory was being brought to its final stage of perfection 
by Laplace, Poisson, Carnot and others, Herapath (1790-1868) in 1820 and 
1821, as well as Waterston (1811-1883) in 1845 proposed their kinetic the- 
ory 19 . However both of them were rejected by the Royal Society. It must be 
emphasized that what needed to be established was not simply a connection 
between heat and molecular motion, for that was already admitted by many 
scientists and was not considered incompatible with the caloric theory; it 
was rather the notion that heat is nothing but molecular motion, and the 
idea that molecules move freely through space in gases rather than simply 
vibrating around fixed positions. This statement could not yet be accepted. 



15 The ratios of volumes of reacting gases are integer numbers. 

16 Let us notice that Avogadro was the first scientist having used the term molecule in 
the modern sense. The word is derived from the french word molecule meaning "extremely 
minute particle" which comes from the Latin word molecula. This term is the diminutive 
of mole meaning "mass, cluster, great number, multitude". On the other hand, the word 
atom is derived from the Latin term atomus used by Lucretius, and from the Greek word 
arofioq, which means uncut. 

17 The atomic weight multiplied by the specific heat of an element is independent of the 
nature of the element. 

18 For further information on the history of atomic theory in chemistry, see for example 
Pullman (1995) and Leicester (1956). 

19 For further informations on their theories, see respectively Brush (1957a, 1957b). 
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During the period between 1842 and 1847 the general scientific and intel- 
lectual climate implied quite simultaneous works by different scientists on a 
concept which is one of the most important reasons of the revival of kinetic 
theory: the conservation of energy. Indeed, the influence of the German 
romanticism and Naturphilosophie, especially the idea that there must be 
a single unifying principle underlying all natural phenomena (Brush, 1967). 
Furthermore, this period inherited different discoveries of various conver- 
sion processes 20 and gave birth to this principle usually attributed to Mayer 
(1814-1878) ( [Mayer, 18421 ) and Joule (1818-1889) d Joule, 1847| . Mayer em- 
phasized the philosophical generality of the principle while Joule provided 
the experimental verification in particular cases. From that moment heat, 
mechanical work, electricity and other apparently different entities are con- 
sidered as different forms of the same thing, now called energy. And in 1847 
Joule and Helmholtz (1821-1894) indicated quite clearly that mechanical en- 
ergy is regarded as the basic entity. It was this prejudice toward mechanical 
explanations that made the kinetic theory appear to be an obvious conse- 
quence of the principle of conservation of energy. If we are convinced that 
heat and mechanical energy are interconvertible, what is more natural than 
to conclude that heat is a mechanical energy? 

The real breakthrough for the kinetic theory took place when Kronig 
(1822-1879) assumed in 1856 (Kronig, 1856 ) that the molecules of gas move 
with constant velocity in straight line until they strike against other molecules, 
or against the surface of the container 21 . But the kinetic theory was still 
confronted to objections like those supported by Buys-Ballot (1817-1890). 
Since the kinetic theory claimed that the velocities of molecules were of the 
order of several hundred meters per second, he pointed out that we would 
expect two gases to mix with each other very rapidly. However, the expe- 
rience shows that the process takes a certain time, of the order of several 
minutes. In order to answer this objection he showed that, in real gases, 
the molecules could not travel without colliding with other molecules. Con- 
sequently, Clausius (1858) introduced a new concept: the so-called mean 
free path of a molecule between two successive collisions. This concept was 
major not only for the further developments, but also because it establishes 
in concrete terms one of the most fundamental ideas of the kinetic theory 



20 Oersted's discovery of electromagnetism (1820), Seebeck's discovery of thermo- 
electricity (1822), Faraday's many discoveries in electricity and magnetism, and many 
others. 

21 Actually his publication did not represent a real advance compared to previous works 
by Bernoulli and Herapath. But his influence in the physics community was important 
and induced an special interest among the physicists for the kinetic theory. 
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of gases rejected in the past: that molecules can move freely through space 
and yet collide with each other. 

The early kinetic theorists assumed that molecules tending to equilib- 
rium move all at the same velocity. Maxwell (1831-1879) was the first scien- 
tist who introduced the idea of random motion for the molecules - hence sta- 
tistical considerations 22 . In 1860, in his first paper on kinetic theory entitled 
Illustrations of the dynamical theory of gases (Maxwell, 1860), he suggested 
that, instead of tending to equalize the velocities of all molecules, the succes- 
sive collisions would produce a statistical distribution of velocities in which 
all might occur, with a known probability. For thermal equilibrium he could 
then derive from symmetry considerations his famous distribution function 
which in modern notation is given v 2 ex.p(—mv 2 /2kBT). 

Hence statistical character appeared to be a fundamental element of kinetic 
theory. 

In 1859, Maxwell came to the kinetic theory as an exercise in mechanics 
involving the motions of systems of particles acting on each other only by im- 
pact. Being interested in viscosity he proposed a mechanism which allowed 
him to establish a relation between the (shear) viscosity and the mean free 
path. By considering gas divided into parallel layers and supposing that the 
motion is uniform in each layer but varying from one to another, he showed 
that the viscosity n should have to be proportional to the mean free path 
(I), the mass density p and the molecular mean velocity (v), relation written 
as 

VM = \{V) p (I) . (11) 



He also published this result in his first article on the kinetic theory ( Maxwell, 1 860). 
However, as Clausius and himself concluded, the mean free path is inversely 
proportional to the density, which implied the surprising result that the vis- 
cosity in a gas does not depend on its density. Whereas it was known that 
this dependence does exist in liquid, and according to the opinion of Stokes, 
Maxwell thought his predictions were absurd and therefore that the kinetic 
theory is wrong, or at least inadequate. But at that time a few experiments 
on the gas viscosity had already been realized. The first studies probably 
began in the 1840's with e.g. Graham (1846, 1849). However, accurate ex- 
periments had not yet been done on the viscosity of gases, therefore, in 1866, 



22 Actually Maxwell was influenced by works on statistics, in particular the works re- 
alized by Quetelet (1796-1874) on the height distribution of a population of soldiers. 
Quetelet stated the idea of the 'average man' with 'average height' etc., as the result of 
an experimentally established normal distribution function. Maxwell worked out the idea 
of an 'average molecule' in analogy with that of the 'average man'. 
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Maxwell himself carried out his own experiment and found out that viscosity 
remains constant over a large range of pressure (Maxwell, 1866). This work 



played an important role in the development of the kinetic theory and its ac- 
ceptance by most of scientists who so far had been doubting it. For instance, 
in 1865, Loschmidt (1821-1895) gave the first convincing estimate of the di- 

o 

ameter of an air molecule (about 10 A which is about four times too large) as 
well as the Avogadro number (about 6.025 x 10 23 ) ( |Loschmidt, 1865 ). On the 



other hand, Maxwell formula, combined with the equation of state for real 
gases derived in 1873 by van der Waals (1837-1923) 23 (|van der Waals, 1873D, 



belongs to the list of thirteen different phenomena mentioned in 1913 by Per- 
rin (1870-1942) and allowing one to evaluate the atomic magnitudes as well 



as the Avogadro number (Perrin, 1991 ). 



As the kinetic theory establishing a microscopic mechanical basis to the 
macroscopic processes such as viscosity was attracting more and more at- 
tention by physicists, the idea to reconcile the second law of thermodynam- 
ics with the principles of mechanics emerged and the first scientist having 
been concerned by it was Boltzmann (1844-1906). His first major achieve- 



ment (Boltzmann, 1868) was to extend Maxwell's distribution law to the 
case of an external force field (field deriving from the potential energy V) 
is present. He obtained the so-called Boltzmann factor which, combined 
with Maxwell's velocity-distribution law, constitutes the basic principle of 
statistical mechanics (Maxwell-Boltzmann distribution law). 

In 1872, Boltzmann attempted to establish an equation describing the 
changes in the distribution resulting from collisions between molecules. In- 
deed he considered the single particle distribution function f(r,v,t) so that 
f(r,v,t)5rSv gives the average number of molecules in the infinitesimal vol- 
ume 5r5v around the position r and the velocity v. The evolution of this 
distribution function is governed by the integro-differential equation which 



is now called the Boltzmann equation ( [Boltzmann, T8721 iBoitzmann, 19951 ) 



Here Jb is the binary collision term 24 taking only two particle collisions 

23 Work based on the virial theorem introduced by Clausius (1870) : "the mean kinetic en- 
ergy of the system is equal to the mean of its virial multiplied by — i; E = — | Y^i r » ' Fj", 
where r ; and Fj are respectively the position of the ith particle and the force acting on it. 

24 Jb(/, f)=J °"dift |vi — v| (f'f{ — ffi) dvidQ. where / and /i, /' and f[ are the func- 
tions of the two particles respectively before and after the collision, crdiff represents differ- 
ential cross section and u> the angular variables of integration. For a detailed derivation 
of the Boltzmann equation, see . Boltzmann (1995) and Dorfman (1999). 
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into account, which is a good approximation for dilute gases. An important 
assumption made by Boltzmann is the so-called Stosszahlansatz 25 , which as- 
sumes that the velocities of colliding particles must be uncorrelated. Later 
Jeans (1877-1946) developed this statement and introduced the assumption 
of molecular chaos 2 ® QJeans, 1903} |Jeans, 1904| ). Boltzmann showed that 
collisions always push / toward the equilibrium Maxwell distribution. In 
particular the quantity H = J f log / always decreases with time unless / is 
the Maxwell distribution in which case H remains constant. This is Boltz- 
mann 's H-theorem < (Bol tzmann, 1872 1. Boltzmann suggested that 
H could be considered as a generalized entropy having a value for any state, 
contrary to the thermodynamic entropy defined only for equilibrium states. 
On the other hand, he explained the reason why the Maxwell distribution 
law is the one corresponding to the thermal equilibrium by showing that 
this distribution is the one most likely to be found, because it corresponds 
to the largest number of microstates. By comparing the two approaches 
(the kinetic approach and the one based on probabilities), Boltzmann could 
conclude that the process of irreversible approach to equilibrium, which is 
a typical example of entropy-increasing process, corresponds to a transition 
from less probable to more probable microstates. Entropy itself can therefore 
be interpreted as a measure of probability. By defining f2 the probability of a 
macrostate, which is proportional to the number of associated microstates, 
Boltzmann proposed a generalized entropy of a macrostate S = kslogO,, 
now called the Boltzmann entropy. 

With the turning of the century Gibbs (1839-1903) transformed the ki- 
netic theory into a more abstract mathematical theory based on the con- 
cept of ensembles of similar systems. Actually Boltzmann (1884) himself was 
working in this way by introducing the so-called Ergoden, which corresponds 
to the microcanonical ensemble of Gibbs (1839-1903). In his definitive trea- 
tise, Gibbs (1902) presented this ensemble as one of a series of ensembles. 
This new formalism appeared to be a powerful technique for calculating and 
predicting the physical properties of bulk matter at equilibrium from the 
properties of atoms. 

Despite all the brillant discoveries due to the atomic hypothesis of matter 
accumulated at the end of the 19th century, a lot of scientists (amongst the 
most important) continued to reject the atomic theory. In general, they be- 
longed to movements called equivalentism represented by especially Berth- 



25 "Assumption about the number of collisions" . 

26 For some informations on the distinctions between both assumptions, see 
jEhrenfest and E hrenfest , 19900 . 
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elot (1827-1907), energetics with Ostwald (1853-1932) and Duhem (1861- 
1916), or empiriocritism, the doctrine of Avenarius (1843-1896) and Mach 
(1838-1916) 27 . The debate then were still vigorous between anti-atomists 
and the partisans of the kinetic theory, especially Boltzmann. Fortunately, 
during the first decade of the 20th century, the triumph of the kinetic theory 
and the atomic hypothesis arrived with the work done in 1905 by Einstein 
(1879-1955) on the Brownian motion (Einstein, 1905). Being aware of the 
attacks on kinetic theory by Mach, Ostwald and others, Einstein started 
his article by contrasting the predictions made by thermodynamics and by 
the kinetic theory (in particular by pointing out that thermodynamics dis- 
tinguishes the Brownian particle and the "hypothetic" molecules compos- 
ing the liquid, whereas kinetic theory does not and says that the colloidal 
particle should follow the Clausius' equipartition theorem). By combining 
Stokes' formula for the force on a sphere moving through a viscous fluid, 
and the formula for the osmotic pressure of dissolved molecules, he derived 
an expression for the mean-squared displacement of the Brownian parti- 
cle submitted to the molecular agitation, which gave a way to evaluate the 
Avogadro's number Na v - In 1908, Perrin (1870 - 1942) made different exper- 
iments in application to the Einstein's Brownian-motion theory and found 
Nav = 6.7 x 10 23 flPerrin, 19081 - Later, in 1913, he published his now famous 
book Les atomes ( |Perrin, 1991 1 in which he proposed thirteen different ex- 
periments based on different phenomena for evaluating Avogadro's number. 
The quite narrow range 6.0 — 7.5 x 10 23 thus obtained is such that we now 
attribute the demonstration of the discrete feature of the matter to Perrin, 
especially due to his efforts to convince the opponents of atoms' existence 28 . 
Whereas Ostwald and others recognized it, some of them like Mach and 
Duhem continued to reject them (Pullman, 1995). 

If the performances made by the kinetic theory were already numerous, 
new works and predictions were added and thus confirmed the depth of this 
theory. In 1911, Enskog (1884-1947) predicted that diffusion would appear 
in a mixture when one imposes a temperature gradient, phenomenon now 
called thermal diffusion. One year later, Chapman (1888-1970) derived an 
expression for the thermal diffusion coefficient. 

Around the end of the 19th century and the beginning of the 20th cen- 
tury different works based on the mean-free-path method were achieved by 

27 For further informations, see for example Pullman (1995), Brush (1967) and Kubbinga 
(1990). 

28 Let us mention that the recent measurements gives us 6.02214199(47) x 10 23 for the 
current value of the Avogadro's number l |Mohr and Taylor, 1998Q . 
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a lot of physicists (Maxwell, Boltzmann, Tait, Rayleigh,...) 29 to improve the 
expression for the viscosity in low-density gases. But one has to wait for the 
almost simultaneous works carried out by Chapman in 1916 and Enskog in 
1917 to obtain an expression for the viscosity derived from the kinetic equa- 
tion 30 ( Chapman, 1917} Enskog, 1917 ) . The derivation known today is that 
of Enskog, and is based on a series solution of Eg, 1)12(1 for the distribution 
function by introducing a parameter A. He therefore derived analytical ex- 
pressions for the transport coefficients for general potential of interactions. 
In the case of the hard-sphere potential, the viscosity coefficient is written 
as 

5 mkRT 
77s = 1.0162— J—^- 13 
16 er z V vr 

where a and m are respectively the diameter and the mass of the particles. 
This expression for the viscosity is now called the Boltzmann viscosity. This 
new expression confirms the temperature-dependence and the absence of 
density-dependence predicted by earlier theory, especially that of Maxwell 
in Eq, (|ll|) 31 . Furthermore, being available for any potential interaction, the 
Chapman-Enskog's developement turned out to be a method for determin- 
ing the actual force law of real molecules by comparing the experimental and 
theoretical results for transport properties as viscosity and thermal conduc- 
tivity. This was the way followed by Lennard-Jones (1894-1954) until the 
early thirties when he proposed the now called Lennard-Jones 6-12 poten- 
tial which is still considered today as the most realistic potential existing 
between particles (Lennard-Jones, 1931 ) 32 . 

Maxwell and others showed both theoretically and experimentally that 
the viscosity coefficient is independent of the density in dilute gases. How- 
ever, even in moderately dense gases, it appears that this property is no 
longer observed. In 1922, Enskog proposed an extension of his previous 
method in order to predict this dependence (Enskog, 1922). Considering a 
hard-sphere system he assumed that the collision rate in a dense gas would 
be changed by a factor x which could be related to the equation of state. He 
modified the Boltzmann equation (|12|) by replacing ffi in which both / and 



29 Especially the dependence in temperature and the factor | in Eg. 1111 . 

30 Notice that Chapman started from the Maxwell's transport theory whereas the En- 
skog's derivation was based on Boltzmann's. For a brief discussion on both methods, see 
Brush (1962a). 

31 By substituting the mean velocity TT and the mean free path 1 by their expression in 
terms of temperature and atomic properties, one obtains r\ = ^ 

mts of the development of expre 

see for example Brush (1970) 



32 For more historical elements of the development of expression for interatomic forces, 
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/i are evaluated in the same point in space, by f(x,y,z, ...)fi(xi, y%, z%, ...), 
where the point points (x,y,z) and (xi,yi,Zi) must be separated by the 
distance a. In such a way, he obtained 33 




) 



(14) 



in which the factor x = 1 + §6o™ + 0.2869 (b n) 2 , b 




and n the num- 



ber density. A density-dependent expression for the viscosity is therefore 
obtained 34 . One had to wait for a decade to compare Enskog's theoretical 
predictions to experimental data. Indeed, Michels and Gibson (1931) mea- 
sured the viscosity of nitrogen gas at pressure up to 1000 atm (1013250 hPa) 
and obtained a good agreement with the semi-empirical Enskog theory pre- 
sented below. Using this method, they also evaluated the dimensions of the 
gas molecule at different temperatures. Let us point out that the only work 
allowing Enskog to test his theory was the one accomplished by Warburg 
and von Babo (1882) for carbon dioxide. The agreement is quite good (see 



In the fifties and sixties, new experimental data for noble gases (for vis- 
cosity as well as for thermal conductivity) were added to allow a comparison 
and a test of theoretical predictions 35 . The particular advantage of Enskog's 
theory compared to other theories for dense gases is that it requires the ad- 
justment of only a few parameters like the diameter of the spheres in order 
to observe an agreement between Enskog's viscosity and the experimental 
data for a larger range of density 36 . 

33 For the development of Enskog's theory, see for example Chapman and Cowling (1970) 
and Hirschfelder et al. (1954). 

34 Let us mention that, in 1899-1900, by considering the analogy with the modification 
of the equation of state for ideal gas when the effect of finite molecular size is taken into 
account, and the correction the mean free path due to the effect of excluded volumes on 
the collision rate, Jager (1865-1938) modified the dilute gas viscosity formula I1H and 
obtained 



where A = 1 + — -p- + ■ ■ ■ . The similarity with the Enskog viscosity 11146 is remarkable 
( |Brush, 19761 - 

35 see for example Senger (1965, 1968), Hanley et al. (1972). 

36 Indeed effective values may be attributed to the parameters b and %■ K was Michels 
and Gibson (1931) who introduced such a procedure by identifying the pressure P in the 
equation of state for rigid spheres with the thermal pressure T(dp/dT)v of the real gas 



Table H. 
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0.719 
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0.704 


0.705 


0.590 


0.722 


0.724 


0.660 


0.756 


0.754 


0.730 


0.795 


0.792 



Table 2: Kinematic viscosity ^ [10 3 poises /{g/ cm?)] versus the mass density 
[g/cm s ] at temperature T = 313.451T, near the critical point. Comparison 
between the experimental results obtained by Warburg and von Babo (1882) 
for carbon dioxide and Enskog's predictions. As the density increases, he 
kinematic viscosity goes through a minimum after which it increases. En- 
skog's predictions confirm this property (data given in Enskog (1922)'s pa- 
per). 



In Tabled it is shown that the bulk viscosity is vanishing in the dilute 
monoatomic gases (argon and helium) as predicted by the Enskog theory. 
Indeed Enskog's bulk viscosity (|15|) for hard spheres 

C = 1.002x^(6 n) 2 (15) 

vanishes for very low density. However, a non- vanishing ratio C/v nas 
been measured in liquid argon fNaugle and Squire, 1965 ) as well as in dense 



gaseous argon (Ma digosky,T967| ). A comparison with the modified Enskog 



theory briefly presented above can be found in Hanley and Cohen (1976). 
Consequently Stokes' relation is not shown to be justified against many 
cases, not only through theoretical predictions, but also experimentally. 



For low densities, one should require that the the Enskog viscosity 114H reduces 
to the Chapman-Enskog viscosity 1131 by requiring lim n ^o X ~ 1, so that we get 
( Hirschfcld er et al., 1954) 

b = B + T— 
dT 

where B is the second virial coefficient. 
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In the sixties it was observed that transport coefficients, in particular 
the viscosity coefficient rj, cannot be expressed in a power series in terms 
of the density n (77 = t\b + f]\n + r]2n 2 + r/3n 3 + . . . ) Indeed it has been 
shown that correlations between molecules is observed over large distances, 
larger than the range of the intermolecular interaction. Hence Dorfman and 
Cohen (1967) obtained theoretically that the coefficient of the quadratic 
term contains a contribution proportional to the logarithm of the density so 
that the previous density expansion should be rewritten as 

rj = rjB +7/1 n + r/ 2 ti 2 In n + 7/3 ra 3 + . . . . (16) 

Nevertheless although the importance of the discovery of long-distance cor- 
relations is great in the modern kinetic theories, their contributions do not 
seem to be important, as it was showed by Sengers (1966) for the viscos- 
ity in a hard-disk gas , and later by Kamgar-Parsi and Sengers (1983) in a 
gas of hard spheres . On the other hand, since the early seventies several 
attempts have been made to detect this logarithmic-density dependence of 
experimental data, in particular for the viscosity coefficient 37 . Even if the 
viscosity is generally the transport coefficient measured with the highest 
precision, it has never been shown that the addition of the logarithmic term 
was necessary. Although we now know that Enskog viscosity is not exactly 
correct, the latter may still be used in order to compare numerical results. 
This is why, despite these facts, we consider Enskog's theory in this thesis 38 . 

Following Enskog theory proposed in 1922, we have to wait 1946 to 
observe a revival of interest for the extension of kinetic theory to dense 
gases and liquids. Indeed, after an original idea of Yvon (1935), a hierarchy 
of coupled equations ruling the time evolution of the reduced n-particle 
distribution functions was derived by Bogoliubov, Born, Green, Kirkwood 
and Yvon such that, by truncating this BBGKY hierarchy, kinetic equations 
were derived leading to expansions of the transport properties in terms of 
the particle density ( Yvon, 1935} Bogoliubov, 1946] Born and Green, 1946} 



Kirkwood, 19 46). 



In 1928, Pauli (1900-1958) derived a master equation in order to apply 



quantum mechanics to irreversible processes (Pauli, 1928). This derivation 
also contained an assumption analogue to Boltzmann's Stosszahlansatz ~ the 
repeated random phase assumption. Since then, different works have been 
completed by means of master equations, such as the generalization of the 



37 See for example Kestin et al. (1971, 1980, 1981). 

38 A historical survey of the discovery of the logarithmic-density dependence is given by 
Brush (1972), pp.72-80. 
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classical Boltzmann equation to quantum systems by Uehling and Uhlenbeck 
(1933). In the fifties, master equations were derived for weakly coupled 
systems by van Hove (1955, 1957, 1959) as well as Brout and Prigogine 
(1956). Kinetic equations were also obtained for plasmas since thanks to 
Landau (1936) and Vlassov (1938). In this context, Balescu (1960) and 
Lenard (1960) derived in the sixties the so-called Lenard- Balescu transport 
equation for plasmas. 

In parallel, the relationship between molecular fluctuations in systems 
due to thermal agitation, and transport coefficients characterizing the irre- 
versible processes was studied starting with Einstein's work on the Brownian 
motion (Einstein, 1905), giving the self-diffusion of the Brownian particle in 



terms of the variance of its position 

x(t) - x(0)f 

D = lim — — '- , (17) 

t— »oo 2 t 

to reach the linear response theory developed in the fifties by Green (1951, 
1960), Kubo (1957), and Mori (1958a). This theory expresses the transport 
coefficients in terms of the integral of autocorrelation function. In particular, 
the viscosity coefficient is written as the autocorrelation of the microscopic 

Current J{t) = £ fl ^PaxPay + \ Ea E&^a F x ab) (r aby) 

r) = lim —3— [°° (J(O)J(i)) dt . (18) 
y^oo Vk B T J 

Moreover, aiming at establishing an Einstein-like relation for transport 
processes, Helfand (1960) derived formulas relating the transport coefficients 
to the variance of their associated quantities now called Helfand moment 
G(t). In particular, in the case of viscosity, such relationship is written as 

\G(t)-G(0)] 2 ) " 
rj = lim lim x 1 ; G(t) = 2> aa; (t)ya(t) , (19) 

a=l 

and has the same structure than Eq. (|17|) . On the other hand, it can be 
shown that the last two relations are equivalent by noting that J(t) = d< ^u ■ 
The transport coefficients are usually calculated based on some kinetic 
equation. Each kinetic equation involves a stochastic assumption such as 
Boltzmann's molecular chaos which is reached by some truncation of the 
evolution equation. Such phenomena as the transport processes, playing 
such an important role in nature, must present a more fundamental justi- 
fication in terms of the intrinsic properties of the underlying microscopic 
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dynamics. We can conclude that microscopic chaos seems to be the key for 
a better understanding of the irreversibility of the macroscopic phenomena. 



5 Theory of chaos 

It is interesting to note that the idea of considering a gas as a disor- 
dered state is old. During the 16th century Paracelsus (1493-1541) studied 
gases and observed that the air and vapors did not have a fixed volume 
or shape. He then gave them the name of chaos in analogy with x ao( > m 
the greek mythology, being at the origin of the universe. Later, van Hel- 
mont (1577-1644), who is the first to have made the distinction between the 
different gases and the air, invented the word gas according to the flemish 



pronounciation of chaos (Leicester, 1956) 



In the end of the 19th century was generally accepted that a large 
number of degrees of freedom is a necessary condition to observe unpre- 
dictable behavior. This great number was the necessary element in order 
to introduce a statistical description of phenomena having a determinis- 
tic underlying dynamics. However, in 1892, Poincare (1854-1912) showed 
that a low-dimensional deterministic system, namely the three-body prob- 
lem, exhibits an unpredictable behavior which is now commonly called chaos 



( Poincare, 1892 ). Since his work, a large number of nonlinear dynamical sys- 
tems, the time evolution of their variables x being given by -j± = Fi(xj), 
have been discovered. In these systems, chaos appears to be a generic phe- 
nomenon rather than an exotic one (Nicolis, 1995|). 



As a matter of fact, in 1873, Maxwell himself had already emphasized 
the sense of unstable behaviors 

When the state of things is such that an infinitely small variation of the 
present state will alter only by an infinitely small quantity the state at some 
future time, the condition of the system, whether at rest or in motion, is 
said to be stable; but when an infinitely small variation in the present state 
may bring about a finite difference in the state of the system in a finite time, 
the condition of the system is said to be unstable, ft is manifest that the 
existence of unstable conditions renders impossible the prediction of future 
events, if our knowledge of the present state is only approximate, and not 



On the other hand, in the late years of the same century, Hadamard (1865- 
1863) studied the geodesic flows on a negative-curve surface and concluded 
that 



39 Cited by Hunt and Yorke (1993). The text of the conference of Maxwell may be found 
in Garnett (1882). 
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each change of the initial direction, as small as possible, of a geodesic which 
remains at a finite distance is enough for implying an arbitrary variation of 
the behavior of the trajectory ( Hadamard, 1898). 

By considering Hadamard's work in a book, Duhem (1906) also concluded 
in the chapter entitled Exemple de deduction mathematique a tout jamais 
inutilisable that, in Hadamard's billiard, the trajectory obtained mathemat- 
ically becomes unusable for a physicist since the an experimental measure 
realized by any physical procedure as precise as possible is always endowed 
of an error which grows with time. Furthermore, Poincare clearly explained 
that 

A very small cause which escapes our notice determines a considerable effect 
that we cannot fail to see, and then we say that the effect is due to chance 
jPomcare, 1908) , 

implying that determinism and chance may be combined thanks to the un- 
predictability. The only mathematician having taken into account what we 
now call the sensitivity to the initial conditions, property put in evidence 
by Poincare and the others, was Birkhoff (1884-1944) who developed in the 
twenties the theory of billiards ( Birkhoff, 1927| ). Otherwise, one had to wait 



for the works done in particular by Lorenz (1963) in meteorology to con- 
sider again the sensitivity to the initial conditions defining chaos, and the 
attractor he put in evidence. Others like Smale (1967) considered global 
behavior of phase space rather than a particular trajectory, and invented 
his horseshoe (see Fig^) in order to get a visual analogy to the sensitive 
dependence to the initial conditions. Finally, Feigenbaum (1978) discovered 
the universality of the period- doubling bifurcation cascade 40 . 

In 1892 Lyapunov (1857-1918) defended his doctoral thesis The general 
problem of the stability of motion ( |Lyapunov, 1892 - ] ). In his work he proposed 
a method which provides ways of determining the stability of sets of ordinary 
differential equations. He introduced the well-known Lyapunov exponents, 
quantities characterizing the growth (A.; > 0) or decay (Aj < 0) rate of the 
distance between two nearby trajectories in the phase space of dynamical 



systems (Eckmann and Ruelle, 1985). The positivity of such quantities ex- 
presses a dynamical instability and induces a sensitivity to initial conditions, 
property defining chaos. Two trajectories initially very close separate ex- 
ponentially in time. This sensitivity to initial conditions limits the possible 
predictions on the trajectories because they are only known through a given 
precision e^ia!- By considering a rate of separation of very close trajecto- 
ries given by the maximum Lyapunov exponent A max , the error between the 

40 For further informations see Gleick (1987) and Ott (1993). 
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predicted and the actual trajectories grows as et ^ einitial exp(A max i). After 
a finite time, the error becomes larger than the final allowed precision efi na i, 
which then defines the Lyapunov time t Ly ap - (1/A max ) In (efinaiAinitiai)- 
Given these initial and final precisions, predictions after the Lyapunov time 
are no longer relevant. This result is the requirement of going a statistical 
description. 

Inspired by the ideas of Krylov (1979), who argued that trajectories in 
phase space for simple fluid separate exponentially, Sinai and his coworkers 
proved the ergodic hypothesis for billiard systems, and characterized the 
stochastic-like behavior of these deterministic systems (Sinai, 19701. This 
type of stochastic-like behavior was then studied extensively for simple sys- 
tems with few degrees of freedom. 

The unstable character of dynamical systems is therefore such that even 
deterministic systems can generate random behaviors. But on the other 
hand, this instability of the dynamics produces information in time allowing 
to reconstruct the system trajectory in phase space. Indeed the separation 
in time of nearby trajectories gives us the possibility to distinguish the tra- 
jectories. In this context, in the late fifties Kolmogorov and Sinai applied 
the concept of entropy per unit time introduced a decade before in the in- 
formation theory by Shannon to the dynamical systems. He defined the 
so-called Kolmogorov-Sinai entropy /iks (Eckmann and Ruelle, 1985). This 



bend this way 



Figure 1: Smale's horseshoe: Succession of stretchings and foldings illus- 
trating the evolution of phase-space due to the chaotic dynamics. 
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new quantity measures the (exponential) rate at which information is ob- 
tained in time in random processes . The KS entropy per unit time is the 
dynamical analogue of the entropy per unit volume defined in equilibrium 
statistical mechanics. Later, in the seventies, Sinai, Bowen, and Ruelle pro- 
posed the mathematical foundations of a new formalism for chaotic systems 
applying the techniques used in statistical thermodynamics, which is called 
the thermodynamic formalism ( Ruelle, 1978| Beck and Schlogl, 1993 ). In 



this formalism, a new quantity, the topological pressure, was defined and 
played a role in dynamical systems very similar to that of the free energy 
for statistical-mechanical systems. 

The two aspects of chaos, that is the dynamical instability characterized 
by the Lyapunov exponents, and the dynamical randomness by the KS en- 
tropy are strongly related, the second one being the consequence of the first 
one. In 1977, Pesin proposed a theorem for closed systems (systems with- 
out any escape of trajectories out of the initial phase space) which is now 
known under the name of Pesin's identity, which relates the KS entropy to 
the sum of all the positive Lyapunov exponents of the system ( Pesin, 1977} 



Eckmann and Ruelle, 1985 ) 



h K s = J2 Xi ' (20) 

A,>0 

As we saw, the instability of the dynamics in phase space induces an 
exponential separation of trajectories characterized by the Lyapunov expo- 
nents. Because the phase-space volume accessible to the trajectories is finite 
(e.g. the volume defined by the energy of the system), they have to fold onto 
themselves. We then observe successive stretchings and foldings of phase- 
space volumes, as the Smale's horseshoe illustrates it (see FigQ), creating 
strange objects (after an infinite number of such operations) nowadays called 
fractals. 

Whereas the term fractal is associated with Mandelbrot, the history of 
this intriguing discipline of mathematics began during the 1870's, when con- 
tinuous functions without derivatives were discovered. For a long time, the 
idealization of nature implied a smooth and regular representation of real 
objects (Chabert, 1994). In mathematics, one dealt with continuous func- 



tions such that a tangent could be drawn at (almost) each point. However, 
Riemann (1826-1866) already claimed a contradictory opinion. And in the 
1870's Weierstrass (1815-1897) gave an example having no derivative in any 
point 41 . In geometry, mathematicians, such as Koch (1870-1924) in 1904, 

41 The paper was read in 1872 in the Royal Prussian Academy of Sciences, but was only 
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Figure 2: Koch's curve: the first four steps of construction. The operation 
consists in dividing each segment into 3 parts, and on the one in the medium 
is constructed an equilateral triangle. 



proposed continuous curves without a tangent at any point obtained by an 



elementary geometric construction (Koch, 1904). Koch's curve depicted in 



Fig|2] is a clear example of such objects. Moreover, the Japanese mathe- 
matician O. Takagi (1875-1960) working at Gottingen proposed in 1903 a 
simple example of continuous but nondifferentiable function known today as 



Takagi 's function (Takagi, 1903). On the other hand, studying the Brow- 
nian motion, Perrin observed experimentally that the trajectory drawn by 
the Brownian particle is highly irregular and he deduced that 

If the functions with derivative are the simplest ones, the simplest to be 
treated, functions with derivatives are the exception; or, to use geometric 
language, curves with no tangent at any point become the rule (from Ref. 
( jPerrin, 1991| ). 

Indeed according to the precision of the measure, the length of the trajectory 
is different, the larger the precision is, the larger is the trajectory. The rele- 
vance of the notion of length thus vanishes and has to be replaced by a new 
quantity characterizing such objects, that is the dimension. In this context, 
Hausdorff (1868-1942) proposed in 1919 a new notion of dimension which is 



no longer a whole number, but can take noninteger values (Ha usdorff, 1 919). 
For example, in the case of Koch's curve, its Hausdorff dimension is neither 
equal to one, that is the dimension of a line, nor equal to two, the dimension 
of the plane, but between both of them, that is 1.26 . 

By introducing the term fractal 42 , Mandelbrot (1975) brought back into 
fashion all the old works on these strange objects. He used them in a lot of 
different disciplines such as in physics (turbulence), geography ("How long 

published on the original version in 1895 flWeierstrass, 1895) . 

42 This word is derived from the Latin word fractus meaning irregular, broken. 
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is the Coast of Britain?"), astronomy (craters of the moon), biology, etc. , 
in order to show that objects without a tangent at any point are nothing 
but the rule in nature. The property being common to the different objects 
he treated, and thus the one implying a certain unity in nature, is the self- 
similarity. A fractal object is such that an enlarged part is similar to the 
whole object, which implies a scale invariance . 



6 Irreversibility and microscopic chaos 

In the context of statistical mechanics establishing relationships between 
the microscopic and macroscopic scales, we deal with microscopic chaos im- 
plied by the defocusing character of successive collisions between atoms and 
molecules. Let us consider the Brownian motion of a colloidal particle. The 
high-dimensional microscopic chaos of the surrounding fluid, characterized 
by the spectrum of Lyapunov exponents, induces a dynamical randomness 
given by the KS entropy, which is calculated by Pesin's identity (T20H . This 
huge dynamical randomness appears to be at the origin of the erratic motion 
of the Brownian particle. Consequently it gives us a new interpretation 44 
of the observed stochastic processes in terms of a high-dimensional chaos in 
the microscopic Newtonian dynamics. Therefore the statistical mechanics of 
irreversible processes no longer needs stochastic models in order to describe 
the macroscopic irreversibility. The microscopic chaos provides determin- 
istic systems which present a stochastic-like behavior without need of any 
stochastic assumption. Moreover, chaos and the sensitivity to initial condi- 
tions implying unpredictability for long times justifies the use of statistical 
mechanics, even for systems with a low-dimensional phase space. 

In this perspective, during the last two decades, different theories have 
been developed establishing connections between irreversible phenomena 
(transport and chemical-reaction processes) and the chaotic properties of 
the microscopic dynamics. Among them we find the thermo stated- system 
approach developed in the eighties by Nose (1984a, 1984b), Hoover (1985, 
1991), Evans and Morriss flEvans et al., 1990[ |Evans and Morriss, 1990D , 



and others. 

In this approach, nonequilibrium systems are defined as systems com- 
posed of particles submitted to interparticle forces, to external forces, but 



43 Several fundamental papers on the development of theory of fractals are reprinted in 
Elgar (1993). 

44 An experimental work has been done exhibiting the chaotic character of the micro- 
scopic dynamics by Gaspard et al. (1998). 
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also to a fictitious nonHamiltonian force modeling the coupling to some 
hypothetical thermostat (Evans and Morriss, 1990). For instance, in order 
to study viscosity, the idea is to reproduce a Couette flow induced by a 
shearing force (see Fig. |3J). Hence a velocity gradient is established in the 






t3 dVx 







Figure 3: Illustration of the shearing of the system inducing a velocity gra- 
dient 

dy 

system. However, this gradient leads to considerable viscous heating of the 
fluid and the energy of the system does not remain constant. To deal with 
this problem, it is necessary to introduce an internal thermostat - a ficti- 
tious frictional force - in order to keep constant the energy. Formally, for a 
3-D system, in addition to the shearing force, we have to introduce a term 



with a thermostating multiplier a ( Evans and Morriss, 1990 ) 



dqi Pi , 

-tt = — + IVi 
at m 

dpi j 

dt 



jPyi - api , (21) 



where 7 is the vector (-^,0), and v is the mean velocity. However, such 
dynamical systems violate Liouville theorem which asserts that the phase- 
space volumes are presented by the microscopic dynamics. This violation 
leads to fundamental problems for defining an entropy in nonequilibrium 
steady states. It thus appears as an artefact of a nonHamiltonian force 
that the phase-space volume visited by the system decreases in time and is 
expressed by the non-zero sum of Lyapunov exponents Yli=i < 0. This 
phase-space contraction is introduced in the system through the presence of 
the thermostat (and the quantity a) which takes away the energy given to 
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the system by shearing 7. It can be shown that the relation between the 
viscosity and the maximum and minimum Lyapunov exponents becomes for 
large systems flEvans et al., 1990) 



r](N, 7) = 3nk 2 BT [A max ( 7 ) + A min ( 7 )] , (22) 

7 

the N dependence of rj and A disappearing. The shear viscosity occurring 
in the Navier-Stokes equation is given by rj = lim 7 _>o 77(7). This expression 
relates the Lyapunov exponents to the viscosity coefficient because of the 
violation of Liouville's theorem by the artifitial nonHamiltonian systems. 
Therefore, this method cannot be used for Hamiltonian systems. Several 
works have shown that other ways exist to maintain a system out of equi- 
librium for instance by stochastic boundary conditions or by deterministic 
scattering (Klages~et al., 2000) in which cases the relation ((2*2*|) do not ap- 



ply. Moreover, violating Louville's theorem creates problems in defining the 
entropy for nonequilibrium steady states. 

Another theory developed in the nineties is the so-called escape-rate for- 
malism (Gas pard, 1998} Dorfman, 199 9) proposed in 1990 by Gaspard and 



Nicolis (1990) for diffusion. As for the previous method, the escape-rate 
formalism introduces nonequilibrium conditions. However, instead of modi- 
fying the Hamiltonian equations, we open the system in order to generate an 
escape process. More precisely, we impose absorbing boundary conditions at 
the statistical level of description, keeping the Hamiltonian character of the 
equations of motion themselves. The so-called escape rate is related to the 
studied transport coefficient on the one hand, and to the chaotic quantities 
of the microscopic dynamics on the other hand. 

Let us illustrate the escape process by taking a simple example: the 
diffusion of a particle in the periodic Lorentz gas (see Fig. 0] (a)). This 
system consists of a particle of nonvanishing mass diffusing thanks to the 
successive elastic collisions with a set of two-dimensional fixed hard disks 
forming a hexagonal lattice. The escape process is introduced by limiting 
the scatterers in a direction, e.g. in the x-direction. It is said that the 
particle escapes when \x\ > h where L is an arbitrary parameter. The 
number of particles lying between the boundaries exponentially decreases in 
time at a rate given by the escape rate (see Figs. HJ. 

In 1995, Dorfman and Gaspard extented the escape-rate formalism to the 
other transport processes ( Dorfman and Gaspard, 1995} Gaspard and Dorfman, 1 995) 



Instead of imposing an external constraint like a shearing, as for the pre- 
vious method, the nonequilibrium conditions are introduced by imposing 
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L 




Figure 4: Lorentz gas with absorbing boundary conditions separated by 
the distance L. (a) Escape process of some particles after diffusion in the 
scatterer. (b) Exponential decrease in time of the number of particles lying 
into the limits defined by — ^ < x < 



absorbing conditions 

-f < G xy {t) < +| (23) 

characterized by an arbitrary parameter Xi over the evolution the Helfand 
moment G xy (t) written in Eq. H19j) . When the Helfand moment reaches 
these absorbing boundaries, the system is said to be escaped. In conse- 
quence, considering a statistical ensemble of systems, an escape process 
appears, which is characterized by the escape rate. Since the system is 
in nonequilibrium, the support of the invariant measure is no longer the 
whole phase space, but a fractal repeller 45 . This fractal is constructed with 
the trajectories for which the Helfand moment remain forever into the ab- 
sorbing boundaries (|23[). The escape rate plays the role of intermediate 
between the viscosity as a transport coefficient, and quantities of the mi- 
croscopic chaos. Indeed, on the one hand, the viscosity can be written in 

45 The fractal object is a repeller because of the instability of the dynamics due to its 
chaotic character. 
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terms of a phenomenological escape rate by considering the macroscopic es- 
cape process. On the other hand, the theory of dynamical systems gives a 
relationship between the escape rate and quantities such as the Lyapunov 
exponents and fractal dimensions, typical quantities of chaos ( |Ruelle, 1978} 
Eckmann and Ruelle, 1985} Beck and Schlogl, 1993 1 . By identifying both 



escape rates, a relationship between the viscosity and the microscopic chaos 

is established and is written as ( Dorfman and Gaspard, 1995} Gaspard and Dorfman, 1995 1 



r, = lim (*) 52(1- di) A, , (24) 

Ai>0 

with the Lyapunov exponents Aj and the associated partial dimensions di 
of the fractal repeller. Hence we can show how the microscopic chaos, 
responsible of the fractal character of the support of the invariant mea- 
sure, controls the viscosity ( |Viscardy and Gaspard, 200 3), or also the diffu- 



sion fGaspard and Baras, 1995] ) as well as the reaction-diffusion processes 
( Glaus and Gaspard, 200Tj|Glaus et al., 2004 ). 



The third approach is the hydrodynamic-mode method ( |Gaspard, 1998D 
developed during the nineties by Gaspard and coworkers ( |Gaspard, 1993} 



Gaspard, 1996). As discussed in Section 12 hydrodynamics describes the 
macroscopic dynamics of fluids in terms of equations governing the evolution 
of mass density, fluid velocity, and temperature, such as the Navier-Stokes 
equations (JIJ and the diffusion equation 

^ = DV 2 n . (25) 
dt ' 

Thanks to the kinetic equation developed by Boltzmann nonequilibrium 
statistical mechanics is able to derive these phenomenological equations, 
while the Boltzmann equation is itself derived from Liouvillian dynamics, 
using the Stosszahlantsatz. The solutions of Ea.(|25jl called hydrodynamic 
modes are of the form 

nk(r, t) = exp(ski) exp(zk • r) , (26) 

each mode being characterized by a wavenumber k. The hydrodynamic 
modes are spatially periodic of wavelength A = 2n/k with k = ||k||. They 
decay exponentially in time because the corresponding eigenvalues are real 
and negative 

Sk = -Dk 2 . (27) 

The paradox between the reversible dynamics at the molecular level and 
the irreversible macroscopic phenomena that was exposed in section[2]can be 
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studied in this context. Indeed the exponential decay of the modes 1)26(1 to- 
ward the thermodynamic equilibrium seems apparently incompatible with 
the Hamiltonian dynamics governing the motion of atoms and molecules, 
which is reversible and preserves phase-space volumes. Therefore, until 
recently the hydrodynamic modes were not described in terms of the Li- 
ouvillian dynamics, but at the intermediate level of the kinetic equations. 
Here, we briefly present a recent method based on the construction of hy- 
drodynamic modes of diffusion in terms of the microscopic deterministic 
dynamics. 

Let us consider the deterministic Liouvillian equation for Hamiltonian 
systems 

& -l = {H, f} = Lf, (28) 

and the solution of this equation is 

/t(X) = F f /o(X) (29) 

where P t = exp(Li) is the so-called Frobenius- Perron operator. 

Boundary conditions are required to solve the Liouville equation. This 
method does not use the simple periodic boundary conditions. Although 
we consider here iV-particle systems periodically extended in position space 
and forming a lattice £ N , such as the Lorentz gas, the probability density is 
allowed to extend nonperiodically over the whole lattice so that the periodic 
boundary conditions does not apply and have to be replaced by the so-called 
quasiperiodic boundary conditions 46 . A Fourier transform must be carried 
out in position space to reduce the dynamics to the cell at the origin of the 
lattice (1 = 0). A wavenumber k is introduced which varies continuously in 
a Brillouin zone reciprocal to the lattice 

A,k(r, p) = Y, ex PH k • !)/*( r + 1. P) • (30) 
ieC N 

In particular the hydrodynamic mode of wavenumber k is an eigenstate 
of the operator T\ of translation by the lattice vector 1 

Ti* k = exp(ik • l)¥ k . (31) 

The wavenumber k characterizes the spatial periodicity of the observ- 
ables and of the probability densities. Each Fourier component of the proba- 
bility density evolves differently in time, which requires the introduction of a 

46 If the time evolution of the coordinates F = (r , p) = (ri , rjv, pi , pjv) of the N 
particles is governed by a first-order equations T = F(T), the vector field F is therefore 
symmetric under discrete position translations: F(r, p) = F(r + 1, p) with 1 G £ N . 
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new Frobenius-Perron operator .R k depending explicitly on the wavenumber 
k. As the operator of translation T\ commutes with the Frobenius-Perron 
operator i? k we can find an eigenstate common to both the spatial transla- 
tions and the time evolution 

# k tf k = exp(s k i)¥ k . (32) 

At vanishing wavenumber, we recover the dynamics with periodic bound- 
ary conditions which admits an invariant probability measure describing the 
microcanonical equilibrium state. In contrast, an invariant probability mea- 
sure no longer exists as soon as the wavenumber is non-vanishing. Instead, 
we find a complex measure which decays at a rate given by the so-called 



Pollicott-Ruelle resonance s k (Gaspard, 1998 ). This measure defines the 



hydrodynamic mode of wavenumber k and the associated Pollicott-Ruelle 
resonance s k gives the dispersion relation of the hydrodynamic mode l|27|). 

Such microscopic hydrodynamic modes present an important difference 
with the phenomenological hydrodynamics. This difference holds in the 
fact that they are mathematical distributions or singular measures. The 
impossibility of constructing eigenstates in terms of functions has its origin 
in the pointlike character of the deterministic dynamics and in the property 
of dynamical instability. Therefore, in order to have a representation of the 
eigenstates we have to consider its cumulative function 

F k (0 = f * k (r c 0< (33) 

Jo 

where T^r is a curve of parameter ^ in the phase space. 

The singular property of the eigenstates plays a fundamental role in the 
understanding of irreversible processes. Indeed it leads to the result that the 
cumulative function -F k (£) is fractal. In the case of diffusion in the periodic 
Lorentz gas, thanks to the thermodynamic formalism, it has been shown 
that the diffusion coefficient D is related to the Hausdorff dimension du of 
the cumulative function i*k(£) ( |Gaspard et al., 2001 ) as 



D = A lim ^SOzl (34) 
k^o k 2 

where A is the positive Lyapunov exponent. Let us mention that a similar 
study has been done in reactive-diffusion systems (Claus and Gaspard, 2004). 



Such processes like diffusion and viscosity are typical irreversible phe- 
nomena. This irreversibility can be expressed by the production of entropy. 
It has been shown that the fractal structure of the diffusive modes plays 



37 



a crucial role in the positivity of the entropy production in Hamiltonian 
systems (Ga spard, 1997{ Gilbert et al., 2000 ). Indeed when the distribution 
functions are smooth, there is no change in the Gibbs entropy and no posi- 
tive irreversible entropy production. The presence of the singular character 
therefore appears to be the fundamental element for an understanding of 
the second law of thermodynamics in terms of fractals. 



7 Conclusions 

In this paper, we have outlined the history of the origin of nonequilibrium 
statistical mechanics, treating the paradox between the reversible character 
of the microscopic dynamics and the irreversible transport processes, such 
as the viscosity. After having reviewed the development of hydrodynam- 
ics and the concept of viscosity, we have shown the development of the 
kinetic theory of gases mainly since Clausius and Maxwell who aimed to 
describe macroscopic phenomena in terms of the microscopic dynamics, and 
how appeared the first interpretation of the origin of irreversibility given by 
Boltzmann. Boltzmann's arguments were based on the probability that a 
system be in a certain microscopic state. For instance, consider an ideal gas 
in a box. Since the state in which all the particles are distributed in the 
whole vessel is infinitely more probable than the one in which the particles 
are concentrated only in the half part, it therefore seems natural that the 
system must tend to the homogeneous state and that we cannot observe the 
reverse case. Consequently, this interpretation is mainly probabilist. 

We have presented an alternative description which has been developed 
in the last few decades. The rise of the theory of chaos and its application in 
dynamical system theory provide a dynamical intepretation of irreversibility 
in terms of the underlying microscopic chaos. Indeed, the typical instability 
of the dynamical chaos induces randomness which is a main characteristic 
of the transport processes (think of the erratic movement of a Brownian 
particle in a fluid). The successful results accumulated since the eighties 
confirm that microscopic chaos is an essential element at the origin of the 
irreversibility. 
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